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Abstract 

The kinetics of a periodically driven nonlinear oscillator, bistable in a nearly resonant 
field, has been investigated theoretically and through analogue experiments. An activa- 
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tion dependence of the probabilities of fluctuational transitions between the coexisting 

i-C ■ attractors has been observed, and the activation energies of the transitions have been 

O ■ 

calculated and measured for a wide range of parameters. The position of the kinetic 
• i— i . 

phase transition (KPT), at which the populations of the attractors are equal, has been 
established. A range of critical phenomena is shown to arise in the vicinity of the KPT 
including, in particular, the appearance of a supernarrow peak in the spectral density of 
the fluctuations, and the occurrence of high-frequency stochastic resonance (HFSR). The 
experimental measurements of the transition probabilities, the KPT line, the multipeaked 
spectral densities, the strength of the supernarrow spectral peak, and of the HFSR effect 
are shown to be in good agreement with the theoretical predictions. 
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I. INTRODUCTION 

The physical kinetics of bistable systems presents a number of intriguing phenomena 
related, in particular, to the existence of substantially differing characteristic relaxation 
times. One (or several) of these, t r , characterise (s) the "thermalization" of a system in 
the vicinity of stable states, i.e. local relaxation and fluctuations about these states. The 
others characterise ffuctuational transitions between the stable states. They are given by 
the reciprocal transition probabilities W^ 1 (the indices i, j enumerate the stable states), 
and usually 

Wijtr < 1 (1) 

The concept of bistability is meaningful, provided only that (1) is fulfilled, in which case 
a system will spend most of the time fluctuating about one or other of the stable states. 
If the parameters of the system pass through the range of bistability in a time much less 
that W^ 1 , the system will display hysteresis: it will tend to remain within one or other of 
the stable states, depending on the prior history. For fixed system parameters, however, 
over times ~ W^ 1 a stationary distribution over the stable states is built up and the 
system "forgets" which of the stable states was occupied initially. 

For thermal equilibrium systems the transition probabilities are usually given by the 
Arrhenius law, W oc exp(— E a /T), where Tis the temperature and E a is the characteristic 
activation energy of the transition. In the case of a Brownian particle the quantity E a is 
simply the depth of the potential well from which the particle escapes [1]. For nonequi- 
librium systems the calculation of the transition probabilities is a nontrivial problem. A 
rather general approach to its solution has been proposed for dynamical systems driven 
by external Gaussian noise (see [2] for a review). In this case W oc exp(— R/a) where a 
is the noise intensity, while R is given by the solution of a certain variational problem. 

In the general case of a bistable system, the characteristic activation energies Ri and 
R 2 for the transitions 1 — > 2 and 2 — > 1 differ one from another. Consequently, for 
sufficiently weak noise, i.e. for small a (when (1) is fulfilled and in addition IF is much 
smaller than the reciprocal correlation time of the noise), the probabilities W 12 and W 2 i 
of the transitions 1 — > 2 and 2 — ► 1 differ exponentially. The latter is the case also for the 
stationary populations w\ and w 2 of the stable states, 
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For most parameter values of the system, the ratio w\jw 2 is either exponentially small (for 
R 2 — Ri 3> a) or large (for R x — R 2 ^> a) and the system occupies with an overwhelming 
probability the states 2 or 1 respectively. Only within a very narrow range of parameters 
where \R\ — R 2 \ ~ a are the populations w\ and w 2 of the same order of magnitude. 
In this range a kinetic phase transition occurs: as we discuss below, the behaviour of a 
noise-driven dynamical system is to some extent analogous to that of a thermodynamic 
system with coexisting phases {e.g. liquid/vapour) within the range of its first order phase 
transition, where both phases are well manifested (with comparable molar volumes, for 
example). 

A well-known signature of systems experiencing phase transitions is the strong asso- 
ciated increase of fluctuations. In the case of dynamical systems, not only will there be 
small fluctuations about the stable states, but there will also be large fluctuations related 
to transitions between the states (the analogue of the fluctuational creation of macrobub- 
bles in a liquid/vapour system). The characteristic time for such fluctuations is obviously 
the reciprocal transition probability W~ x . It is quite natural to expect that these large 
and very slow fluctuations will give rise to intense and extremely narrow (with a width 
~ W oc exp{—R/a)) peaks in the susceptibility of the system and in the spectral density 
of fluctuations (SDF) [3]. For a Brownian particle fluctuating in a symmetric double- well 
potential (that is, exactly at the phase-transition point, w\ = w 2 ), the corresponding 
peak in the SDF at zero frequency has already been observed [4]. The exponentially 
fast broadening of this peak with increasing noise intensity gives rise [5] to low frequency 
stochastic resonance [6-8] i.e. to the increase with increasing noise of the signal and the 
signal/noise ratio in a system driven by a low- frequency periodic force. 

An important class of bistable systems is those displaying bistability when driven by an 
intense periodic field, but which are monostable otherwise. (Note that some such systems 
may also display multistability, and/or dynamical chaos, when subjected to even stronger 
periodic fields.) A variety of them are investigated in, for example, nonlinear optics (in re- 
lation to optical bistability, see [9]). The different stable states here correspond to periodic 
attractors with differing amplitudes, phases (and sometimes frequencies) of constrained 
vibrations. A well-known example [10] of such a system is the nonlinear oscillator driven 



by a nearly resonant field. This system is interesting not only as an archetypal model for 
the investigation of a periodic-field induced bistability, but also because it refers directly 
to a peculiar and interesting physical system, a relativistic electron trapped in a magnetic 
field and excited by a resonant cyclotron radiation [11]. 

In addition to its markedly nonequilibrium character, which provides a good test for 
theories of fluctuational transitions in thermally nonequilibrium systems, the model of 
a resonantly driven nonlinear oscillator also enables one to investigate specific phenom- 
ena arising in systems with coexisting states of forced vibration within the range of the 
kinetic phase transition. Since there is a special frequency in such systems, namely the 
frequency of the external field up, the fluctuational transitions between the stable states 
should modulate the response of the system at frequency u>f' extremely tall and narrow 
spectral peaks near ujp are therefore to be expected, both in the susceptibility [2] and 
in the SDF [12]. Such supernarrow spectral peaks have been observed in an electronic 
analogue experiment [13]. Because the widths of such peaks increase extremely rapidly 
(exponentially) with noise intensity, it is to be anticipated that the signal/noise ratio for 
a signal at a frequency close to ujp will also increase with increasing noise intensity, i.e. 
that there will be a manifestation of high frequency stochastic resonance (HFSR). 

In the present paper we present detailed results of our investigation of the features of 
the SDF in a periodically driven system, including the onset of a supernarrow spectral 
peak in the region of the kinetic phase transition. In Sec II below the theory of kinetic 
phenomena for a periodically driven oscillator is given, including the results of a numer- 
ical calculation of the "activation energies" of the fluctuational transitions and explicit 
expressions for the spectral density of fluctuations and for the generalized susceptibili- 
ties. In Sec III the experimental simulation of the oscillator by an analogue electronic 
circuit is described. In Sec IV the theoretical and experimental results are compared with 
each other, and the new critical phenomena, the onset of the supernarrow spectral peak 
and of the high-frequency stochastic resonance, are discussed. Sec V contains concluding 
remarks. 



4 



II. THEORY OF THE SPECTRAL DENSITY OF FLUCTUATIONS OF 
AN OSCILLATOR BISTABLE IN A PERIODIC FIELD 

A. The Model 

In this section we explore the behaviour of a nonlinear oscillator subject to the com- 
bined influences of a periodic field Fcosupt and a weak random force f(t). The equation 
of motion of the particular oscillator considered (single-well Duffing oscillator) is 

q + 2Tq + uj\q + 7 g 3 = Fcoscu F t + /(*) (3) 
The oscillator is assumed underdamped and the periodic field nearly resonant 

T, \5u\ <C lofi 5uj — ujf — ojq (4) 

A characteristic amplitude of vibration for which the oscillator will obviously be strongly 
nonlinear, the nonlinear length l n , is determined by the condition that the nonlinear term 
7g 3 in (3) should be as large as the linear one, so that l n = (wo/|7|)5. If the amplitude F 
of the periodic force is sufficiently large that vibrations of amplitude a ~ l n are excited, 
then the oscillator (3) in the absence of the random force f(t) is known [14] to display 
deterministic chaotic phenomena (see also [15]), with the boundaries of the domains of 
attraction to various attractors often being fractal [16]. 

In the case of an underdamped oscillator, strong nonlinear effects can also arise for 
much smaller values of F, for which the vibration amplitudes a are correspondingly much 
less than l n [10]. This is because the eigenfrequency of the vibrations of a nonlinear 
oscillator depends on their amplitude, u e ff = cu(a), and it is the interrelation between the 
detuning of the field with respect to the eigenfrequency and the damping, \ujp — u(a)\/T, 
that determines the strength of the response. For small T, the latter ratio can vary 
markedly with a, even while a <C l n , and it is this feature that can give rise to the 
coexistence of different stable solutions for the amplitude a. We may note that, for 
a <^ l n , the vibrational amplitudes at the overtones are ~ o 3 //^ a, and dynamical 
chaotic phenomena do not occur. Because the noise intensities to be considered here are 
relatively weak, the system seldom strays far from the attractors, and practically never 
goes as far as l n ; the probability of finding it there is exponentially small compared to the 
probabilities of transitions. 



5 



Under conditions for which |7|(g 2 ) <C wf, the motion of the oscillator consists of rela- 
tively fast oscillations with slowly varying amplitude and phase. The characteristic scale 
for these variations is determined by the friction coefficient Y and the detuning 5u of the 
field frequency oof with respect to the oscillator eigenfrequency uoq\ the characteristic scale 
for the coordinate q is (7 /ojf8oj\~z . In describing the "slow" motion it is convenient, in 
the spirit of a standard averaging method (cf [17]) to transform to the rotating frame. We 
thus change from q, q to the complex dimensionless envelopes u, u* and the dimensionless 
time r, 

T = \6u\t 

The equation of motion in terms of the variable u following from (3) - (5) takes the form 

^ = v + V f(r) (6) 

v = v(u,u*) = —r)u + iu{\u\ 2 — 1) — i(3^ 

where 



are respectively the reduced damping coefficient and the dimensionless field intensity. The 
equations (7) as written corresponds to the particular case 



5uj > 0, 7 > (8) 

The generalization to the case where the signs of Sou and 7 are arbitrary is straightforward. 
We note that bistability can occur only for ^ySuj > 0; simultaneous change in the signs of 
7 and 5uj will result in mirror reflection of the spectra considered below with respect to 

Up- 
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The random force /(r) appearing in (6) is proportional to /(f) in (3), 



and, if /(%) is Gaussian white noise of characteristic intensity B, such that 



(/(f) /(f)) = 2TB5(t-f) (10) 
then /(t) is asymptotically a two-component white noise, 



(/»/>')> = (/»/>')> = 



</(r)/V)) = 4a*(r - -0 (11) 

Ct = 5"= 

i6w|.r 

The correlator (Re/(r)Im/(r')) is fast-oscillating; the slow variables u,u* therefore per- 
ceive the components Re /(r), Im /(r) as independent white noises of equal intensity; 
which is why the correlators (/(t)/(t')) and (/*(t)/*(t')) are set equal to zero in (11). 
We also note that the relations (11) can be asymptotically fulfilled even where the initial 
noise f(t) is not 5-correlated; it suffices that its correlation time be small as compared 
with the "slow" process times r -1 (but not necessarily as compared with tup 1 [3]). 

The dynamics of the system (6) depends on the values of the three dimensionless 
parameters involved: rj, j3 and a. We shall assume the noise to be weak, so that 

a«l (12) 

To zeroth order in /(r), Eq (6) describes the autonomous motion (note that we consider 
it in the frame rotating with the frequency of the external field) of a system with one 
degree of freedom and, correspondingly, with two dynamical variables u and u* (or Re 
u and Im u). The stationary solutions of the equation du/dr = v give, in accordance 
with (5), the states of steady forced vibration of the oscillator. The values of the complex 
envelope u in the steady states follow from the relation v — 0, and are given by 
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% = V^(N 2 - 1 + ir i)~\ <KM 2 ) = o 

(j>(x) = x(x-l) 2 + rfx-(3 (13) 

where j enumerates the real roots of the cubic equation (13) and can take on the values 
1, 2 or 3. Eq (13) is readily interpreted. As a result of nonlinearity, the frequency of the 
eigenvibrations of the oscillator depends on their amplitude a as 

3 

8 

On substitution of this expression into the well-known expression for the amplitude of the 
forced vibrations of a linear damped oscillator of frequency a; (a) 



[u 2 F - u\a)f + 4c4r 2 
with account taken of the relation a 2 = ^\u\ 2 ujf\^\/\j\, which follows immediately from 
(5), one simply obtains the cubic equation 

0(M 2 ) = o 

In the parameter range where this equation has three real roots, the oscillator is bistable: 
the forced vibrations with the smallest (j = 1 in (3)) and the largest (j = 2) amplitudes 
are stable; there is also the unstable steady state, j = 3, with an intermediate value of 
\uj\ 2 oc a 2 . The phase of the stable forced vibrations is determined by the argument of Uj 
in (13). The range of f3, rj for which (13) has three solutions, and thus bistability occurs, 

(1 2) 

corresponds to the approximately triangular region bounded by the full curves f3 B ' (rj) 
of Fig 1, i.e. the bifurcation curves. 

Thus, as the amplitude of the periodic force is gradually increased from a small value at 
fixed frequency (see e.g. vertical line a-a' in Fig 1), the system moves from monostability 
(one small limit cycle), to bistability (two possible limit cycles of different amplitude), 
and then back again to monostability (one large limit cycle). Some analogy can be drawn 
between the bistability and the liquid/gas coexistence region of a Van der Waals system. 
As the spinode point, which corresponds to the Van der Waals critical point, is approached, 
the two stable (and one unstable) solutions of (13) coalesce and the amplitudes of the 
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large and small limit cycles (liquid and gaseous phases) correspondingly approach each 
other, to become indistinguishable at the spinode (critical point) itself. Consequently 
just as in the Van der Waals case, it is possible to move quasistatically from an initial 
state that is a small limit cycle (gaseous state) to a final state that is a large limit cycle 
(liquid state) without undergoing the analogue of a first order phase transition or passing 
through a mixed-phase coexistence region: all that is necessary is to take a route through 
the parameter space that passes outside the spinode (the critical point). 

The above analysis makes sense provided that the basins of attraction are smooth and 
regular throughout the region of phase space likely to be visited by the system. The basins 
have been computed within the bistable regime (see Fig 1) in the absence of noise, and are 
shown as Poincare sections (values of q, q for t = 2-nnuOp 1 +0 O ) i n Fig 2. We emphasize that 
the data in Fig 2 refer to the initial oscillator described by (3) with f(t) = 0. In addition 
to the dimensionless parameters 77, f3, this system is characterised by the parameter Y/uj f 
which, in the present case, was set as Y/up = 0.0184. The results were obtained by the 
usual "grid of starts" method [15], allowing the system to evolve from different starting 
points in the (q, q) phase space and noting in each case the attractor to which it was 
drawn. Thus, all starts in white areas lead to the large amplitude attractor (the • in the 
white area), and all starts in black areas lead to the small amplitude attractor (the • in 
the black area). Fig 2 shows the evolution of the basins with increasing j3 for fixed i] as 
Poincare sections for the same phase. It is intuitively reasonable that the black basin (for 
the small amplitude attractor) should be dominant at small (3, just within the region of 
bistability. As (5 increases, the white basin (for the large amplitude attractor) grows until 
the central regions of the two basins have become equal in area. With further increase of 
(3 the black basin continues to shrink, finally disappearing at the upper boundary (Fig 1) 
of the bistable region. 

The most important feature of Fig 2 for present purposes is that it confirms the basins 
to be (within the range and resolution of the computations) simple, smooth, and regular 
as already stated above. The shapes and positions of the attractors are close to those 
given by the approximate equations of motion (6) in the absence of noise; we note that the 
latter equations do not display chaos or fractal boundaries. Because we are interested in 
the regime of weak noise intensity for which the system only makes occasional transitions 
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between the attractors, we can be confident that it spends almost all of its time in the close 
vicinity of either one or other of them, and that the probability of fluctuations carrying 
it out to regions of phase space where the basins might be irregular or fractal (far beyond 
the range plotted in Fig 2) is exponentially small. 

Finally in this section, we draw attention to the importance of a slightly different 
model, closely related to (3), that is likely to be more easily realised in experiments on 
systems excited by laser radiation with a randomly varying amplitude: 

q + 2Tq + uj\q + 7 g 3 — [F + f(t)} cos u F t (14) 

The significance of (14) arises because of the very high value of the driving frequencies 
uj f in optical experiments, which means that external noise f(t) introduced from a con- 
ventional noise generator will in practice be far from white; indeed the cut-off frequency 
of the generator is likely to be much smaller that ujf- 

Nevertheless, transforming (14) to the rotating frame again gives Eq (6), except that 
the noise is now given by 

plus a term which varies as f(t)exp(2iujt). The limited spectral width of f(t) that we 
have assumed implies that this term will only have a very small effect on the equation of 
motion of the slow variable u(t). It is the high-frequency components of the noise that 
determine the random dynamics of a nearly-resonantly-driven underdamped nonlinear 
oscillator; the components of the noise with frequencies far from u are filtered out. In 
contrast to /(r) (9), the new /(r) has correlator 

(/»/>')> = (/»/>')> = -°S(t - r') 

(f(r)r(r')} = a5(r-r') 

i.e. instead of two independent components, the new /(r) has only one. Nonetheless, the 
analysis presented below can be easily extended to give similar results for the system (14). 
Modulation of the periodic driving force by noise has pushed the effect of low frequency 
noise into the high frequency range. 
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B. Transition probabilities and the spectral density of fluctuations 

The most obvious effects of noise on the behaviour of the oscillator are, first, the 
onset of fluctuations about the stable states and, secondly, the occurrence of fluctuation- 
induced transitions between the states. Provided that the noise is weak, in accordance 
with (12), the system will spend most of its time in the close vicinity of one of the stable 
states: only very rarely will a sufficiently large fluctuation occur to cause a transition to 
the other stable state. The dependences of the probabilities of the transitions on the 
characteristic noise intensity are of the activation type. 



The activation energy Ri for the transition from state i is given by the solution of a 
variational problem: the corresponding variational equations and the algorithm for their 
numerical solution are discussed in the Appendix. The resultant dependences of Ri on 
P for the lower (i = 1) and higher (i = 2) amplitudes of the forced vibrations in the 
limit of small reduced damping r\ were considered in [3]. Numerical results for four 
values of rj are shown by the circles in Fig 3(a)-(d). It is evident that Ri decreases, 
and i?2 increases, monotonically with increasing (3 i.e. with the characteristic resonant 
field intensity. For the values of (3 corresponding to the upper and lower bifurcation lines 
(3 B 1,2 ^ in Fig 1, Ri and R 2 respectively vanish (as the states 1 and 2 coalesce with the 
saddle point and then disappear). The dependence of Ri on (3 for (3 close to (3b\Q) * s 
universal, Ri = Gi(rj)\f3 — /^(f?)] 3 / 2 , and is shown by the full lines. (The explicit form 
of Gi has been considered earlier, cf [3]). The numerical and asymptotic results are in 
good agreement for not too small rj, where the optimal path of the escape (in the rotating 
frame (5): see Appendix) is not a small-step spiral. For small 77, however, the numerical 
algorithm is not accurate enough and results in the discrepancies seen in Fig 3(a); as 
discussed in the Appendix, the data in this range can better be obtained in a different 
way. The dependence of R12 on the frequency detuning 77 for (3 lying in the central part 
of the interval (Pb\t)), Pb\v)) is rather sharp, especially at small 77 where [3] R\^ oc 77. 
As 77 approaches its critical value (the spinode point in Fig 1) given by 



Wij = const, x exp(— Ri/a) 



(15) 
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the values of decrease rapidly [3], as (77 — i]k)^ 2 - Here, too, the numerical and analytic 
results are in good agreement. It is evident that, as rj approaches t]k, the range of (3 over 
which Ri is well-described by the asymptotic law \(3 — Pb\v)\ 3 ^ 2 increases relative to the 
total range of bistability ~ Pb\v)\- Qualitatively, this is related to the fact that 

the optimal path of the escape is approaching a straight line on the plane (u f ,u"), and it 
becomes nearly straight for all /3, 77 close to the spinode point. 

The analytic and calculated values of R 1 and R 2 allow us to plot on Fig 1 the dashed 
curve (3o(f]) specifying the points at which the activation energies are equal, closely ap- 
proximating the line of the kinetic phase transition (KPT) at which the populations are 
equal, 

R 1 (f3,r ] ) = ifcOM, P = Po(v) (16) 

(Note that the criteria R\ = R2 differs from w± = w 2 only by virtue of variations in the 
prefactor ~ a in (15); the criteria become identical as a — > 0). For parameter values 
far from this curve, the transition probabilities W 12 and W 21 are seen from (15) to differ 
exponentially strongly. Correspondingly, the stationary populations Wi,w 2 of the states 
as given by (2) are also exponentially different, 

Wi/w 2 oc exp[(i?! — R 2 )/a\ 

Only in the close vicinity of the fia{rj) curve will the transition probabilities and stationary 
populations be comparable and it is here, therefore, that one may expect to observe the 
characteristic steady state fluctuation phenomena associated with transitions between the 
attractors; we do not consider here the transient fluctuation effects that arise when, for 
example, the parameters are swept through the bifurcation lines in Fig 1. 

A revealing characteristic property of a fluctuating system is its spectral density of 
fluctuations (SDF). The SDF of the co-ordinate of a periodically driven oscillator, Q(u), 
is given by 

1 roc 

Q{uj) = -Re / dtexp(iut)Q(t) 

7T JO 

Q(t) = r Km ^ f_ T drW + r) - (q(t + r))][g(r) - (q(r))} (17) 
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We note that a periodically driven system is in general nonergodic, so that Q(t) is not 
equal to the time correlation function 

({q(t + r) - (q(t + r))][q(r) - (q(r))]) 

defined in terms of ensemble averaging (....); in fact, the latter quantity oscillates with r 
at frequency cu F , as can be seen from (5), (6); Q(t) actually corresponds to this quantity 
smoothed over r. 

The ensemble-averaged value of the coordinate, (q(t)), is equal to the value of q(t) 
averaged over equal instants of time modulo 2ti/u f , 

N-l 

(q(t)) = lim N' 1 V q(t + 2nnujp 1 ) 

In the case of weak noise, two principal contributions to Q{uj) can be identified [12, 
18]. The first of these arises from small fluctuations about the stable states. It is equal 
to the sum over the states j of the corresponding partial SDFs, Qj(cu), weighted by the 
state populations Wj given by (2) (cf also [19]). The second contribution Qt r {w) comes 
from the (relatively infrequent) fluctuational transitions between the states. Thus (cf [5]) 

Q(u) = Y,WjQj{u) + Q tr (uj) (18) 
j 

To calculate the partial SDF, Qj(cu), for the state j when the noise intensity a is 
small, it suffices: to linearize the v term on the right hand side of (6) with respect to 
small deviations in (u — Uj), (u* — u*)] to substitute the solution of the resultant linear 
equations into (17), taking due account of (5); and to perform averaging. The result is of 
the form 



Au F aT 2 (u - u F ) 2 + 2{uj - c^iy^M 2 - 1) + Y 2 {v) + 2?T>/) 
Q ^ ' ~ ~3\^r [(u - uo F ) 2 - TV. 2 ] 2 + AT\u-uo F ) 2 



v) = 1 + r 2 (3|^| 2 -l)(| Mj | 2 -l) (19) 

where \uj\ 2 for j — 1, 2 is given by Eq (13). It is evident from (19) that Qj(ou) is peaked 
near the frequency uj f of the external field. Its intensity will be proportional to the noise 
intensity. The shape of Qj(uj) will be discussed in Sec IV. 
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The second term in (18), Q tr {oj), can readily be calculated if one notes that the 
populations Wj of the stable states fluctuate in time with a characteristic relaxation time 
(W u + W 2 i) _ \ so that 

= -(w 12 + W 21 ) Wl {t) + W 21 

at 

w 2 (t) = l-w^t) (20) 

(The values of Wj appearing in (2), (18) correspond to the stationary solutions of (20)). 
In the case of weak noise, these fluctuations can be shown to result in a contribution to 
Q(u) of 

2U F \U F -U \ 2 W 12 + W 21 

Qtr{u) = ^-rn \{ui) - {u 2 }\ Wl w 2 — — 2 21 

3vr|7| (W 12 + W 21 ) +{u- uj f ) 2 

Here (u)j denotes the ensemble average value of u in the state j. In the zero noise limit 
(u)j is simply Uj. For the purposes of comparison with experiments performed at finite 
noise intensity, (u)j can be expanded as a perturbation series in the small parameter a. 
To first order we obtain 



(5 Uj ) = {i^ 2 + 3|w,| 4 - Q\u 3 \ 2 + 2] - (3hf - 2)( v 2 + 2\ Uj \ A - 3\ Uj \ 2 + 1)} 

We note that the spectral peak Q tr (uj) is extremely narrow: its width is determined 
by the transition probabilities, so that it is exponentially small and much smaller than 
the damping parameter Y which determines the "dynamical" relaxation of the oscillator 
towards either of its stable states. The product wiw 2 , which determines the intensity of 
Q tr (ou), can be seen from (2), (15) to be exponentially small for almost all values of (3,7], 
with the exception of those within the very narrow range (the phase transition region) 
where wi ~ w 2 ~ 1. Thus, the onset of the fluctuational transition-induced spectral peak 
Qtr{w) is a specific phase transition phenomenon (see Sec IV). 

C. The susceptibility and high frequency stochastic resonance 
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The effect of a weak trial periodic force on thermal equilibrium systems, is the onset of 
vibrations at the frequency of the force; their amplitude is characterised by a susceptibility 
which can be expressed in terms of the SDF via the fluctuation dissipation theorem [20]. 
If the system is being driven by a strong periodic force Fcosu^i, so that it is far from 
thermal equilibrium, the additional weak force Aexp(-iQt) gives rise to vibrations not 
only at its own frequency Q, but also at combination frequencies \u ± u F \, \Q ± 2up\, ■■■ 

We shall consider the linear response of the bistable oscillator to a nearly resonant 
trial force with a frequency Q close to uo,uf'- 

\Q — LJ \, \Q — U)p\ <C CJp 

In this case, a strong response is to be expected. It will be most pronounced at the 
frequency Q and at the nearest resonant combination frequency, which is 2ujf — £1. Thus, 
one can seek the trial force-induced modification of the ensemble-averaged coordinate q 
in the form 



5(q(t)) ~ x(tyAexp(-int) + X (VL) Aex^[i{2u F - Q)t] 



(22) 



That is, we may suppose that the linear response is characterised by two coefficients 
(generalised susceptibilities), x(^) an d X(Q). The absorption/amplification of the trial 
field is characterised by Im x(ty- ft was shown in [3] that, in the vicinity of the KPT, 
interesting features occur in the absorption spectrum. 

To calculate the susceptibilities we transform to the slow variables u, u* (5) in the 
equation of motion (3) with the additional force Aexp(-iflt) added to the right hand 
side. The resultant equations for u, u* take the form 



^ = v(u,u*) + rjf(r) 



du* 



V*(U,U*) + T)f*(T) + lA(r) 



(23) 



Mr) = 
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It is evident from (23) how the second term in (22) arises: it is due to the addition 
oc exp[— i(Q — u F )t] to u, which is then multiplied by exp(icu F t) when q(t) is calculated 
in (5). 

If the random force f(r) is weak, the main effects of the additional term oc A(t) in 
(23) are: (i) to cause small amplitude periodic vibrations of u, u* about their stable values 
Uj, Uj] and (ii) via the change in the probabilities of fluctuational transitions, to modulate 
periodically the populations of the stable states. These effects give rise to expressions for 
the generalised susceptibilities of the form 



x(n) = ^^(Q) + Xtr(n) (24) 
j 

where Xji Xj are the partial susceptibilities related to the corresponding vibrations about 
the stable states, and Xtr, X tr are related to the trial force-induced redistribution over 
the states. 

The partial susceptibilities can readily be calculated by linearising (23) near the stable 
states, yielding 

, . i T — i(Q — up) — i(2\uj\ 2 — l)(up — u ) 

x (Q] . . u fa F z ^ 

j[ } ~ 2u F T 2 v]-2tT(n-Lu F )-(n-Lu F y [ ' 

The effective modulation of the transition probabilities by the trial field A, which gives 
rise to the second term on the right hand side of each of Eq (24), arises when its frequency 
Q is very close to cu F , so that 



\Q — uj f \ T, \uj f — uoq\ 

In this case, the trial field smoothly raises and lowers the effective "barrier" between the 
stable states with the period 2k/\VL — u F \, so that the activation energies -Ri,i?2 of the 
fluctuational transitions vary periodically in time [3]. The corresponding additions to Rj 
are given in the Appendix. In turn, they give rise to periodic additions to the transition 
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probabilities Wij (15) and hence to the population Wj of the stable states (20). The final 
expression for the redistribution-induced additions to the generalised susceptibilities is 



Xtr ^ = v~r, -t(( m )i - ( M ) 2 )( — - — ) 



I i(Q — uj f ) 1 

W12 + W 2 l 



Xtr (Q) = — r^Xtrity (26) 

It is evident from (26) that the susceptibilities Xtr(fy, X tr (£l) are large only within the 
range of parameters (3, r], close to the kinetic phase transition, where the populations 
u>i, w 2 of the stable states are of the same order of magnitude. The characteristic range of 
the frequency Q of the trial field within which these susceptibilities are large is determined 
by the transition probabilities. Consequently, it increases exponentially with increasing 
noise intensity (cf [12]). This property gives rise to stochastic resonance [8] i.e. to an 
increase of the signal/noise ratio (SNR) with increasing noise [21] which, as shown below, 
occurs in the present system for a high frequency signal 11 ~ » T. To calculate the 
ratio, we note from (22) that the signal induced by a real field Acosfit is given by 

5(q(t)) = ARe{ X (tyexp(-int) + X(fi)exp [i{2u F - Q)t]} (27) 

Such a signal corresponds to the appearance of <5-shaped spikes in the power spectrum of 
the oscillator at frequencies Q and (2up — fi). This can be seen from (17) if (27) is added 
to q(t), q(t + r) but not to (q(t)). (The latter quantity is included in (17) to subtract the 
(^-function in Q(Q) at the frequency ujf of the strong field). It is evident from (16), (17) 
that the ratios P and V of the strengths (areas) of the spikes at frequency f2 of the trial 
field, and at the combined frequency {2uj f — ft), to the power spectrum in the absence of 
noise are given by 



S 



Q(ny Q(2u F -n) 



s= 1 -A"\ x (n)\\ s= 1 -A"\x(n)\ 2 



(28) 



17 



It follows from (25) that the "partial" susceptibilities Xj(^)> Xj(Q) are independent 
of noise for weak noise, whereas the partial contributions to the SDF Qj(£l) increase 
linearly with the noise intensity. Far from the phase transition region, therefore, where 
the fluctuational transition contributions to the susceptibilities and SDF are small, the 
quantities P and V decrease with increasing noise. Within the phase transition range, on 
the other hand, for small \Q — oof\ ~ Wij, the main contribution to x(f2), X(Q) and Q(Q) 
comes just from the transitions (26), (21). This is because their ratio to the corresponding 
partial contributions is inversely proportional to the (small) noise intensity parameter a. 
If only Xtr, X tr and Q tr are taken into account in (28), one obtains the corresponding 
quantities P tr , V tr 



According to (29), the quantities P tr ,Vt r are independent of frequency. At the same 
time, they can be seen from (15) to increase exponentially with increasing noise intensity. 
This implies the onset of high frequency stochastic resonance within the phase transition 
range, not only at the frequency of the trial field, but also at the combinational frequency 
(2uf — ft). In fact, the ratio V is rather different from the quantity usually considered 
in the context of SR, because no force is being applied at the frequency \Q — 2u>f\' the 
signal is induced by mixing, in a nonlinear system, of the forces at frequencies fl and ujp- 
In relation to nonlinear optics [22] , the phenomenon can be regarded as a type of highly 
selective, resonant, four- wave mixing (actually, multiwave mixing because the effect is not 
just proportional to the squared amplitude F of the strong field). 

We would emphasize that stochastic resonance occurs only within the phase transition 
region. When the parameters (3, r\ of the oscillator are far from this region, the contribu- 
tions Xtr-, X tr , Q tr to the susceptibilities and the SDF in the absence of the trial force are 
exponentially small: P and V differ markedly from P tr ,Vt r , therefore, and decrease with 
increasing a. The dependences of P and V on a, as given by (18), (19), (21), (24)-(26) 
and (28), will be compared with the results of the analogue electronic experiments in 
Section IV below. 




(29) 



III. ANALOGUE ELECTRONIC EXPERIMENTS ON 



THE PERIODICALLY DRIVEN OSCILLATOR 
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In order to test the theoretical predictions of the preceding section, and to find out 
whether they were applicable to a real physical system described by the model equation 
(3), a series of analogue electronic experiments was undertaken. The basis of the ana- 
logue technique has been described in detail elsewhere [23], together with a discussion 
of its advantages and disadvantages. In essence, it is extremely simple. An electronic 
model of the stochastic differential equation under study is built using standard analogue 
components (operational amplifiers, multipliers etc). This is then driven by stochastic 
and periodic forces, as appropriate, and its response is analysed with the aid of a digital 
data processor. 

The circuit used to model (3) is shown in (slightly simplified) block form in Fig 4. It 
was designed and scaled in the standard [23] way so as to optimise use of the dynamic 
range of the active components. Thus, the actual equation simulated (see Fig 4) was the 
integral form of 

R\ „ . R\ R\ x 3 
RiC 1 R 4 C 2 x + —R4C2X = -—x - — — 

R R 
+-^-F' sin u' F t + — i A' sin n't + r)(t) 
Ri Rq 

with 

R 1 = R A = R 5 = 2R 6 = ^ = ^ = 100A;fi 
R 3 = 2.5MQ 
d = C 2 = InF 
t = R\C\ = R4C2 
2T = R 1 /R 3 
Thus the equation actually simulated was 

x 3 F' A 

t 2 x + 2Ttx = —x 1 sin uLt H sin Q't + rt(t) 

10 10 F 10 /w 

which, after the scaling 

t^t'/r, uj f ^uj' f t, Q^Q'r, F-> if/10, A — > A'/IO 
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goes over into (3) with u = 1, a = 0.1. Provision was made for measuring either the 
coordinate q(t) or the energy E = \q 2 + \q 2 + \^iq 4 . 

The circuit model was driven with a sinusoidal periodic force from a Hewlett-Packard 
3325B frequency synthesizer. Its response, a time- varying voltage representing q(t), was 
digitized (12-bit precision) typically in IK or 2K blocks and analysed using a Nicolet 
LAB80 data-processor; for the experiments on high frequency stochastic resonance, where 
larger q(t) data blocks were required (see below) a Nicolet 1280 data-processor was used. 

As expected, the model was found to display bistability within a certain range of 
forcing amplitude and frequency: its response q(t) in the absence of noise for a given set 
of parameters, shown in Fig 5(a) and (b), could have either of two distinct amplitudes, 
corresponding to the two coexisting periodic attractors discussed in Section II. An inherent 
experimental difficulty of the measurements lay in the accurate determination of f3, on 
account of the \uip — c^o| 3 term in the denominator. Because, for the region of interest, 
the forcing frequency up is very close to the natural (zero amplitude) frequency ujq of the 
oscillator, a very small error (typically ± 1%) in ojq inevitably results in a much larger 
error (typically ± 40%) in the value of f3. For this reason, rather than attempting to 
determine o;o from the nominal component values or, directly, by a resonance experiment, 
its value was established precisely by measurement of the range of bistability at a single 
value of 77. Once this had been done, the rest of the region of bistability could be mapped 
out over the whole range of 77, resulting in the square data points of Fig 6. They are seen 
to be in satisfactory agreement with the theoretical prediction of Section II (full curves). 

The energy E = \q 2 + \q 2 + j^q 4 of the oscillator (apart from the coupling energy to 
the force), measured as a function of q in the absence of noise for each of the attractors, is 
shown in Fig 5(b). Both E and q are periodic functions of time, so that the energy can in 
principle take on several values for any given q, depending on how many times q becomes 
zero during one period 2ti/uf- For the Duffing oscillator in the range of parameters 
considered here, where the nonlinearity is relatively small (see Section IIA), q was zero 
twice within a period and therefore E could take on not more than two values for a given 
q. Since, for small enough nonlinearity, the steady vibrations q(t) have components at 
the odd overtones only, i.e. at the frequencies u>f, 3u>f, 5uf ■ ■ ., there is an additional 
symmetry: q(t + ft/up) = —q(t), E(t + t/Vf) = E(t). Consequently, the energy is a 
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unique function of the coordinate q on an attractor. We emphasize, however, that E is 
not conserved: the oscillator acquires energy from the periodic driving force and dissipates 
it through friction. It is evident from Fig 5(b) that the curvature of E(q) is relatively- 
small, providing a clear indication that the amplitudes of the harmonics of E(t) are also 
small; in the approximation (6), they have been ignored. The smearing (thickness) of the 
E(q) lines in Fig 5(b) is an experimental artefact: the values of E and q were recorded 
at discrete intervals of time, and neighbouring pairs of values have been connected by 
straight lines (cf Fig 7(b) below). 

When noise was applied to the driven oscillator, fluctuations about the attractors 
and occasional transitions between them were observed. Fig 7(a) shows an experimental 
example of one such transition; the corresponding E(q) plot in Fig 7(b) shows the fluc- 
tuations in energy about each of the attractors, yielding an envelope that illustrates very 
clearly the shape of the potential. 

Measurement of the transition probabilities between the attractors was not completely 
straightforward because there was a small region of overlap between them in terms of any 
single variable whether measured, for example, in terms of q(t) or of E(q). Thus the 
determination of sojourn times [24] on either side of a fixed value of q or E would not 
have provided the information sought. Instead, the mean first passage time (MFPT) was 
measured between two pre-set criterion levels in energy, which were outside the overlap 
region and unambiguously within each of the attractors. Fig 8(a) illustrates how the 
apparent MFPT varied when one criterion level was kept fixed within the lower energy 
attractor, and the other level was moved through different values. There is clearly a 
plateau region around the noise-free energy level of the attractor for which the MFPT 
was independent of level setting: all of the MFPT measurements to be reported below 
refer to this region. 

To determine SDFs, a standard fast Fourier transform (FFT) routine was used to 
compute the power spectral density of the fluctuations q(t) — (q(t)). In practice, the 
ensemble- averaged signal (q(t)) was determined in a preliminary experiment for each set 
of parameters, averaging a large number (typically 1000) of blocks of q(t) in order to 
obtain good statistical quality. This was possible because of (q(t)) being strictly periodic, 
with (q(t)) = (q(t + 2it/u>f)), and because the phase of (q(t)) was determined with 
respect to that of the field Fcosu;_pt. The resultant (q(t)) was then subtracted from each 
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subsequent realisation of q(t) before the FFT was applied to find the SDF, which was 
itself ensemble- averaged to produce the final result. 

The experiments on high frequency stochastic resonance involved the application of 
an additional weak trial force Acos(Qt) to the system, with Q very close to the main 
forcing frequency u>f- In order to resolve the expected (see above) responses at fl and 
at \2u>f — ^| from the supernarrow peak at u>f, it was necessary to use a relatively large 
block size, which in practice was set to 8K or 16K using the Nicolet 1280 data processor. 



We now compare the theoretical predictions of Section II with the results of he ana- 
logue experiments described in Section III. As already noted above, all of the main features 
expected on the basis of the theory have been observed in the simulations: for example, 
the anticipated bistability of the oscillator was observed and its range in terms of f3, rj 
was found (see comparison of data points and theory in Fig 6) to be in agreement with 
the theoretical predictions; and with weak noise applied to the system, fluctuational tran- 
sitions were observed to be taking place between the stable states. We now present a 
more detailed comparison of experiment and theory considering, in turn, escape proba- 
bilities, spectral densities of the fluctuations and high frequency stochastic resonance in 
Subsections A, B and C respectively. 
A. Transition probabilities 

To characterise the transition probabilities, the average lifetimes (TJ) of the states were 
measured (with the mean time (T,) from the initially occupied state i being measured in 
the absence of back-flow as described in Section III, so that 



with (dwi/dt)o = —WijWi and Wi(0) = 1). Some typical measurements of the average 
lifetime (Tj) = W^ 1 , on a log plot as a function of noise intensity a, are shown in Fig 
8(b). The fact that the data fall on straight lines confirms that the escape process is 
of the activation type, as expected on the basis of Eq (15); the characteristic activation 
energies i?j of the transitions for given f3, 77 can be obtained immediately from the slope 
in each case. Some experimental values of Ri, obtained in this way from a large number 
of measurements similar to those of Fig 8, are presented in Fig 9. In good qualitative 
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agreement with the theoretical predictions (full curves), the value of R\ for the transition 
from the lower amplitude attractor decreases monotonically with increasing dimension- 
less field intensity f3, while R 2 for the transition from the higher amplitude attractor 
correspondingly increases. Both Ri and R 2 increase with the increase of the frequency 
detuning parameter 77. Note that the experimental errors here are relatively large, due 
to the problem of measuring 77, discussed above, and to the effect of small changes ~ ± 
0.5% in uo with ambient temperature. 

The values of /3, 77 for which Ri = R 2 lie extremely close to those values for which 
(Ti) = (T 2 ) [because the effect of the prefactor in Eq (15) is relatively weak], defining 
the kinetic phase transition. The phase transition points obtained from the experimental 
data (for (Ti) = (T 2 ), crosses) fall close to the position of the theoretical phase transition 
line (for Ri = R 2 , dashed) in Fig 6. The discrepancy between experiment and theory 
becomes somewhat larger near the spinode point K where the system is very "soft" and 
extremely weak noise intensities are necessary to make the transition region sufficiently 
narrow and the phase transition itself sufficiently sharp. The influence of uncertainties 
in the experimental parameters, and of internal noise in the active circuit components, 
becomes even more important here. Taking due account of all these factors, it may be 
concluded that theory and experiment are in satisfactory agreement. 

B. Spectral density of fluctuations 

Experimental measurements of the SDF in the vicinity of the oscillator eigenfrequency 
(histograms) are presented and compared with theory (full curves) in Fig 10(a)-(c). It 
must be emphasized that the measurements refer to the spectral density of fluctuations 
about an ensemble average; the subtraction of (q(t)) from each realisation q(t) prior to 
computation of Q{uS) ensures that, when the system remains on one particular attractor 
throughout, most of the 5-shaped at oo F (which is the Fourier transform of (q(t))) gets 
removed. This is why there is very little sign of a spectral peak at ujf in Figs 10(a) and 
10(c). In the KPT range, however, where jumping occurs between the attractors, the 
general appearance of the spectrum is entirely different. In fact, the most striking feature 
of the spectrum is the supernarrow peak [15] that rises in the phase transition range, 
where (5 ~ fl (r)), as seen in Fig 10(b). Its width is very much smaller than either of the 
widths of the other peaks, or the experimentally determined damping constant T, or the 
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frequency detuning ujp — oj§ (which are all of the same order of magnitude). For small 
noise intensities a, this width is unresolved by the LAB80 data analysis system i.e. the 
peak lies entirely within one "bin" of the data-processor's memory. It was necessary to 
increase the noise intensity substantially in order to spread the peak over two or three 
bins. 

The dependence of the intensity / of the supernarrow peak on the distance from the 
phase transition line was found to be exponential, as shown in Fig 11. This feature can 
readily be understood in terms of (2), (15) and (21). According to (21), for small noise 
intensities where (ui) ~ Ui, 

1= / duQ tr {UJ) = — Ul -U2 W 1 W 2 (30) 

J-oo 3|7| 
Not too far from, but not too close to, the phase transition line where, on the one hand, 
\P ~ Po(v)\ ^ 1 an d ; 011 the other, the transition probabilities W 12 and W 2 \ differ sub- 
stantially from each other, the coefficient Wiiv 2 in (21) should behave, according to (2), 
(15) as 

Wl w 2 oc exp{-|i*i-i^||/3-A)(77)|/a} (31) 

where R' { = (dRi/dp)p = p ( v y Therefore, the dependence of I on (ft — (3 ) should indeed 
be exponential. In the immediate vicinity of the phase transition, where the exponent on 
the right hand side of (31) is of order unity, this exponential dependence will be smeared. 
Such smearing is clearly seen in the experimental data (squares, Fig 10), which are in 
good agreement with (30) as indicated by the crosses (wi,w 2 having been taken from 
independent measurements of the transition probabilities). Agreement with (30) based 
on Wiw 2 taken from the simple estimate (31), as indicated by the full lines, is also good: 
the values of R' 12 i n this case were taken from the slopes of the experimental R\ j2 {f3) 
measurements and the prefactor in (31) was taken to be \ so as to give the correct 
maximum value of W\W 2 at the phase transition point j3 = (3o(v) m t ne limit a — > 0. 
The cusp-like dependence of the intensity / of the supernarrow peak in the SDF is a 
characteristic feature of the peak, which itself represents a characteristic phase transition 
phenomenon peculiar to bistable systems. 

For (3,T) lying far from the phase transition line flo(v)i the supernarrow peak is not 
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seen, but there remain (histograms of Fig 10) the much broader and less intense peaks 
in the SDF corresponding to fluctuations about the stable states. These correspond to 
the partial spectra of the first term in (18) and are well described (full curves) by (19). 
The characteristic feature is that they each (J — 1, 2) display a twin-peaked structure for 
a certain range of parameters. It can be seen from (19) that such structure should be 
at its most pronounced for the spectrum Qj(£i) when ^> 1: that is, for sufficiently 
large frequency detuning \ujp — u \ ^> T. Under these conditions, the peaks appear at 
{uj f — u)q) ~ ±ri/j, and their half-width ~ Y is small compared to the distance separating 
them. The twin-peaked structure can be understood intuitively in terms of the forced 
vibrations at frequency in a given stable state being modulated by the relatively slow 
(characteristic frequency ~ up — w ) fluctuational vibrations about this state. We notice 
that the intensities of the peaks in a doublet differ markedly (parametrically strongly, in 
the case of the small- amplitude attractor), so that the intensity of the second peak for 
the small-amplitude attractor is fairly small. 

In the range of the kinetic phase transition, the partial spectra Qi(Q), Q2(&) are 
superimposed and the supernarrow peak is also present. Thus there can be up to five 
separate peaks in the spectrum. A multi-peaked structure is clearly evident in the results 
of Fig 12, which were recorded for a larger detuning and a smaller T. A satisfactory 
quantitative description of such a spectrum cannot, however, be arrived at on the basis 
of (19), because it is significantly influenced by higher order terms in the noise intensity, 
i.e. by vibrations at the overtones of Tuj, that were ignored in the derivation of (19). A 
detailed investigation of such higher order effects is currently being planned and will be 
the subject of a future paper. 

C. High-frequency stochastic resonance 

In searching for evidence of the predicted high frequency stochastic resonance (HFSR) 
phenomenon, the circuit parameters were: 2r = 0.0397; u Q = 1.00; 7 = 0.1; ujf — 
1.07200; D, = 1.07097; F = 0.068; A = 0.006. The frequencies of the additional weak trial 
force and the main periodic drive were therefore very close to each other. A typical SDF, 
measured for a 16K digitized time series in the memory of the Nicolet 1280 data processor 
with input noise intensity B = 0.040, is shown in Fig 13. The central maximum is the 
supernarrow SDF peak of Fig 10(b), here with its finite width clearly resolved (note the 
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highly expanded abscissa scale). A delta function spike is evident, not only at Q, but also 
at the mirror-reflected frequency (2uf — ^) just as predicted in Section IIC. 

The signal strengths (integrated intensities) S(a), S(a), determined from measure- 
ments of the delta spikes, are plotted (data points) as functions of the reduced noise 
intensity a oc B in Fig 14(a) and (b). It is immediately apparent that there are well- 
defined maxima in the plots of S(a), S(a), thereby confirming the occurrence of HFSR 
in (3). The signal/noise ratio P, V (i.e. the ratio of S, S to the "background" SDF in 
the absence of the trial force) have also been measured. As shown in Fig 15, although 
the statistical quality of the data is somewhat poorer (owing to the additional error in 
the background SDF), there is no doubt that P,V each fall and rise and fall again with 
increasing noise intensity. The fall in P, V with increasing a at small a is, of course, 
a feature that is familiar from earlier calculations and experiments on conventional SR 
involving static bistable attractors; the signal/noise falls initially because, for very weak 
noise, the inter-attractor transitions are too rare to make significant contributions either 
to the susceptibility or to the SDF, whereas the background SDF in the denominator 
corresponding to fluctuations about the attractors steadily increases with noise intensity. 

The theoretical predictions, based on Eqs (18), (19), (21) (24-26) and (28), are shown 
by the full curves in Figs 14 and 15. The agreement is not perfect but (given the problem 
with the determination of f3: see Section III), it is within the experimental error and may 
be regarded as satisfactory. The onset of the observed rise in S, S, P, V occurs at the value 
a of noise intensity for which the width of the supernarrow peak in the SDF becomes 
comparable with the frequency difference \Q — ujf\ (provided that the latter is not itself 
exponentially small; cf [7b] where the position of the minimum of the SNR versus noise 
intensity has been discussed for stochastic resonance in a system fluctuating in a static 
bistable potential). It is the increasing role of fluctuational transitions that is responsible 
for high frequency stochastic resonance. These results, and the good agreement obtained 
with the theory of Section II, demonstrate that HFSR for periodic attractors may be 
perceived as a linear response phenomenon, in very close analogy to conventional SR for 
a static bistable potential [7]. 

An intuitive understanding of the mechanism of HFSR can be gained by recalling that, 
under the conditions considered here with |Q — uj f | very small, the system responds to the 
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trial force almost adiabatically. In terms of the phase diagram Fig 1, the beat envelope 
of the combined main and trial periodic forces results in a slow vertical oscillation of the 
operating point p. When this is set (see line p' — p — p") to straddle the KPT line, which 
was the case for present investigations, and the noise intensity is in the appropriate range, 
the system will have a tendency to make inter-attractor transitions coherently, once per 
half-cycle of the beat frequency. The net effect is to increase the modulation depth of 
the beat envelope of the response, thereby amplifying its component frequencies Q and 
\Q — 2uf\- 

The magnitude of the signal at Q has been measured as a function of distance, ex- 
pressed in terms of /3, from the KPT. The result is shown in Fig 16. It exhibits a fast 
cusp-like (note the log scale) decrease of S as (3 moves away from its critical value, demon- 
strating that, like the associated supernarrow spectral peak (see Section IV B), HFSR for 
periodic attractors has the character of a critical phenomenon, in agreement with the 
theory of Section II. 

V. CONCLUSION 

Studies of a very simple nonequilibrium bistable system - a nearly resonantly driven 
nonlinear oscillator - have enable us to observe and investigate a number of new phenom- 
ena of rather general applicability and, in doing so, to test a theoretical approach to the 
calculation of transition probabilities in noise-driven nonequilibrium systems. The onsets 
of the supernarrow peak in the spectral density of the fluctuations and the corresponding 
peak in the susceptibility, and the phenomenon of high frequency stochastic resonance, 
can all be viewed as examples of critical kinetic phenomena in periodically driven sys- 
tems. They may be used, not only to investigate the character and properties of kinetic 
phase transitions (as here), but also to obtain tunable noise-induced amplification of the 
signal/noise ratio and extremely narrow-band filtering and detection of high-frequency 
signals. 

Finally, it is interesting to note that many of the phenomena discussed above provide 
illustrations of the creative role often played by noise in nonlinear systems. The occurrence 
of the extremely strong and highly selective four-wave mixing, for example, arises because 
of the noise, and does not occur in the absence of noise; the dependences on noise intensity 
both of this effect, and of the other critical phenomena discussed in the present paper, 
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are exponentially sharp. 
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APPENDIX 

The calculation in [3] of transition probabilities for systems driven by Gaussian noise 
was based on an idea [25] due to Feynman. His suggestion was that there was a direct 
relationship between the probability density of the paths of the noise-driven system and 
the noise itself. This inter-relationship allows us to write immediately, to logarithmic 
accuracy, the probability density functional for the paths of the system, and to set up 
the variational formulation giving the most probable paths for first reaching a given point 
in the phase space of the system and for the transitions between the stable states (see 
also Ref 2). In the white noise case under consideration, the "activation energy" Rj 
characterising the transition of the oscillator from the stable state j to the stable state % 
is given [3] by the following variational problem 

^w>(HfiH ( - 4i) 

v = v(u,u*), u(—oo) = Uj, -u(oo) = u s 

where v is defined by (6), and Uj and u s are the values of the "slow" variable u for the 
initially occupied stable state and for the saddle point respectively. The general analysis 
of large occasional fluctuations in systems driven by white noise was given by Wentzell' 
and Freidlin [26]. 

In obtaining a variational (Euler) equation for the problem (Al), u and u* should be 
varied independently. The resulting equation can be seen to be of the form 

_ 2 A2| M | 2 - 1) - V 2 uu 2 + Jp(2\u\ 2 + u 2 — 1 — iff) = (A2) 
ar l dr v 

V 2 EE U 2 (U,U*) = 1+?T 2 (M 2 - 1)(3M 2 - 1) 

together with the equation for u* complex conjugate to (A2). The corresponding equations 
for u' ee Rew, u" = Imu were written down explicitly in Ref [2]. An analytic solution can 
be obtained [2] in some limiting cases. In particular, in the small rj limit, the equations 
describe the conservative motion of an auxiliary system with two degrees of freedom, its 
coordinates being u' ', u" and its velocities ii' and u". The motion can be considered as the 
planar motion of a particle of unit mass and unit electric charge in an electric potential 
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where 



g(u ,u ) = -(u + u — 1) — uyp 

and a magnetic field 7i = [4(V 2 + w" 2 ) — 2] normal to the plane. The potential U (V, w") 
is shown in Fig 17. It has three maxima of equal height (=0). They correspond respec- 
tively to the stable states and the unstable stationary state of the system. Note that we 
are dealing with the auxiliary system. Thus, not minima, but maxima of the potential 
U (V, u") correspond to the stable states of the initial system (a circumstance typical of 
the instanton formulation that we are using). The problem (Al) amounts to finding a 
path that starts on one of the outer maxima and arrives at the maximum corresponding 
to the saddle point. 

The numerical solution in the general case of arbitrary values of (f3, rj) can be simplified 
by the following procedure (cf [27]). Near the stable state j, when \u — Uj\ 1, Eq (A2) 
can be linearized in u — Uj, u* — u*. The solution can then be sought in the form 

u{t) - Uj = 5>$ s) exp(A$ s) r) (A4) 



for r — > — oo, with 



pW A (s) *j ~ L ) - VI ^°\ u j\ -*mi 1 > ( a<\ 

j ~ j 2uj(2\ Uj \* - 1 - i V ) [ ' 

The resultant characteristic equation gives four values for the increment A 

Af 2) = af 2 \ \f> = -a?\ Af = -af (A6) 

where af' 2 ^ are the roots of the characteristic equation for the motion of the oscillator in 
the vicinity of the state j in the absence of the random force (the latter being described 
by the linearized equation du/dr = v), 
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For the stable state, Re aj ' ' < 0; note that this implies that z/| > 0. 

It is evident from (A4)-(A7) that the coefficients Af 2) (and thus Bf' 2) ) in (A4) should 
be set equal to zero; otherwise, the path u(t) will not approach Uj as r — > — oo. At this 
point we have arrived at two independent parameters in (A4): and . However, 
it is their ratio which determines the direction in which the system will move along 
the extreme path (A2); accordingly, it is this ratio that should be determined from the 
boundary conditions u(0) = u s (Al). We thus obtain an implicit equation for the single 
quantity Af/Af. It can in principle be solved numerically by a shooting method (cf 
[28]). This method works most effectively for small damping, 77 <§C 1, when the optimal 
path u'{t),u"{t) is a spiral. 

Here, we adopt a relaxation method [29]. The differential equations are cast in the 
form of nearest neighbour difference equations. The boundary values of u, u* were chosen 
at Uj, u s . (The results were the same to the adopted accuracy when iy was taken instead 
of u s , i.e. the optimal path for the escape from the state j was supposed to start from 
Uj and to arrive at the other stable position u^-j] moreover, it turned out that the path 
found in this way passed through, or fairly close to, the saddle point). A guess at the 
solution was tried, and then it was successively improved by assuming, at each step, that 
the true solution was close to the current one and linearising about the latter. Even for 
fairly different initial guesses, the same final solution u(t) was obtained, implying that the 
method is reliable. The same results were obtained for different integration times, which 
were always very large compared to the characteristic dimensionless times ~ 1,?7 _1 , in 
(A2). 

Solutions were used to compute the action integral (Al). The results are summarised 
in Fig 17. The optimal path is indicated by the o points, which are separated by equal 
intervals of time. The motion is naturally at its slowest (points closest together) on the 
maxima of the potential; because of the effect of the "magnetic field" H, it does not pass 
exactly along the ridges. The advantage of the relaxation method is that it is fairly fast 
and convenient. On the other hand, for small rj where the optimal path is a small-step 
spiral, it is less reliable than the method based on solving Eqs (A4)-(A7) which was used 
in [28]. 

When, in addition to the strong field F cosu F t, the oscillator is also driven by a weak 
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force Aexp(— iQt), the expressions for Rj will change. The additions to Rj can readily 
be found when \Q — ujp\ <C T, because the characteristic time of the motion along the 
extreme path described by (A2) is (i.e. rf x in dimensionless units of r); thus, when 
\Q — lvf\ ^ r, the weak field is not changed while the system is moving along the path. 
It is evident from (22) that the functional which should be minimised to give Rj in the 
presence of the field A is given by (Al) with v* having been replaced by v* + iA(r). To 
first order in A, the change in Rj is thus of the form 



where the integral giving fij is calculated along the extreme path for A = 0. 

The expression for fij can be substantially simplified if one notices that the activa- 



respectively, where ip is arbitrary: such a replacement corresponds simply to a shift of 
the time origin in (5) by ip/cup, which should not influence stationary characteristics of 
the oscillator such as Rj. By differentiating Rj with respect to ip for ip — > 0, one finds 
immediately that Im/Xj = 0, and it is then easy to see that 






tion energy Rj (Al) is unchanged when fli is replaced by P^exp^iip) in (6) and (Al) 




(A9) 
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FIGURE CAPTIONS 



1. Phase diagram for the system (3) in terms of the reduced parameters (7). Within 
the approximately triangular region enclosed by the full lines, the system is bistable, 
with two possible stable limit cycles of different amplitude and phase relative to the 
periodic driving force. Outside this region, the system is monostable. The dashed 
line represents the calculated position of the kinetic phase transition. The cuts a-a' 
and p'-p-p" are discussed in the text. 

2. The evolution of the basins of attraction for (3), computed for rf = 0.072 in Poincare 
section with the same phase as that of the driving field. The white regions show 
the attracting basin of the large amplitude attractor, and the black regions show 
the basin for the small one; each attractor is indicated by a •. The values of f3, 
left-to-right, from top-to-bottom, were: 0.0709; 0.0811; 0.0913; 0.102; 0.112; 0.122; 
0.132; 0.143; 0.153; and 0.163. 

3. Calculated activation energies Ri for transitions between the coexisting periodic 
attractors of (3), as functions of (5 for: (a) rf = 0.033; (b) rf = 0.100; (c) rf = 
0.200; (d) rf = 0.333. The circle data points were obtained by minimisation of the 
action energy integral Rj in Eq (Al); the curves are derived from analytic expansions 
[3] around the relevant birfurcation points. In each case, the falling data and curve 
represent R\ and the rising data and curve represent R 2 . 

4. Block diagram of the analogue electronic circuit model of Eq (3). 

5. (a) Variations of the coordinate q(t) with time t, measured for the electronic circuit 
model of (3) within the bistable regime in the absence of noise with (3 = 0.0607, rf = 
0.033. Digitized time series for the small and large amplitude attractors are shown 
superimposed; note the phase difference between them. In (b) the instantaneous 
energy E of the system is plotted as a function of q: the lower curve is for the small 
amplitude attractor and the upper curve is for the larger one. 

6. Comparisons of the calculated region of bistability (between the full curves: cf Fig 1) 
with that measured for the electronic circuit model of (3) (square data points), and 
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between the calculated (dashed line) and measured (crosses) kinetic phase transition 
line. 

7. (a) Variation of the coordinate q(t) with time t, measured for the electronic circuit 
model of (3) in the bistable regime of the presence of noise with (3 = 0.0607, rf = 
0.033, showing a transition between the attractors (a rare event), (b) The corre- 
sponding variation of the instantaneous energy E with the coordinate q (cf Fig 5 in 
the absence of noise). 

8. (a) Dependence of the apparent mean first passage time on the position chosen for 
the upper criterion level for crossings: the energy of the larger attractor for q = 
in the absence of noise is shown by the dashed line. In practice, all measurements 
were made in the plateau region, (b) Logarithmic plots of the MFPT between the 
attractors measured for the analogue electronic circuit model of (3) as a function 
of reciprocal noise intensity a -1 , for rf = 0.033: + from small amplitude attractor 
with (3 = 0.0888; □ from large attractor with (3 = 0.047; O from small attractor 
with (3 = 0.0734. The fact that the data lie on straight lines demonstrates the 
"thermal activation" character of the transition mechanism. 

9. Values of the activation energies Ri measured (crosses) as functions of (3 for the 
analogue electronic circuit model of (3) with rj 2 = 0.033: (a) for Ri, and (b) for R 2 . 
The theoretical curves are the same as those of Fig 3(a). 

10. Spectral densities Q(ui) of fluctuations measured (histograms) for the analogue elec- 
tronic circuit model of (3) with r] = 0.219, a = 8.69 x 10~ 2 for: (a) (3 = 0.048; (b) 
(3 = 0.078; (c) (3 = 0.150. The full curve represents the theory. 

11. Variation of the intensity / of the supernarrow spectral peak with distance from the 
kinetic phase transition line, measured as a function of (3 for the analogue electronic 
circuit model of (3) for rj = 0.219, a = 8.69 x 10 -2 . The squares represent direct 
measurements; the crosses are derived from (21), based on measured transition rates. 
The full lines also represent (21), but for \n(wiw 2 ) given by (31) with measured 

12. An example of the kind of multipeaked spectral density of fluctuations Q{oj), mea- 
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sured in the range of the kinetic phase transition for the analogue electronic model 
of (3), with relatively large frequency detuning and small damping. The parameters 
were: r] = 0.055; (5 = 0.0303; T = 0.0073; a = 0.238. Five spectral peaks - the 
supernarrow peak and four peaks of the partial spectra - are clearly resolved. 

13. Spectral density Q{uj) of the fluctuations of (3) as a function of frequency u with 
an additional weak periodic (trial) force Acos{VLt + 0) added on the right hand 
side, measured for the analogue electronic circuit model. The contents of each 
FFT memory address are shown as a separate data point on a highly expanded 
abscissa (unsealed experimental unit). A smooth curve has been drawn through the 
background spectrum (the supernarrow peak, which has its maximum at ojf) as a 
guide to the eye; vertical lines indicate the delta spikes resulting from the trial force. 

14. The intensities S and S of the 5-shaped peaks in the SDF of the analogue electronic 
model of (3) (data points) with (5 = 0.103, r\ = 0.266 induced by a weak trial 
force A cos Qt, plotted as a function of noise intensity a, compared with theory (full 
curves): (a) at the trial force frequency f2; (b) at the mirror-reflected frequency 
(2lu f - Q). 

15. The signal/noise ratios P and V of the responses at /3 = 0.814, rj = 0.236 to the trial 
force at frequencies Q (circle data and associated curve) and (2uf — fi) (squares), 
measured as functions of noise intensity a for the analogue electronic circuit model 
of (3). The curves represent the theory. For noise intensities near those of the 
maxima in P(a) V(at), the asymptotic theory is only qualitative and so the curves 
are shown dotted. 

16. Dependence of the logarithm of the intensity S of the spectral peak induced by the 
trial force for the analogue electronic circuit model of (3) at a fixed noise intensity 
a = 0.076, plotted as a function of (3, varying across the KPT line. 

17. The three-humped potential U (u', u") of Eq (A3) for the auxiliary system, portrayed 
both as a three-dimensional surface and, below it, in the form of a contour plot; the 
contour altitudes are tabulated on the right hand side. (Note, however, that the left 
hand potential maximum is so shallow that it is barely visible in the contour plot). 
The optimal path of the escape from the focus-1 (the small amplitude attractor) 
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is shown by the line of o points. It goes from this focus (the right hand potential 
maximum) to the saddle point (central maximum). On its way from the saddle 
point to the focus-2 (the left-hand potential maximum corresponding to the large 
amplitude attractor) the system moves, with overwhelming probability, along the 
noise-free path. The o points are equally spaced in time, so that the speed with 
which the system is travelling along different elements of the path may be inferred 
from the density of the points; it moves most slowly at the three potential maxima. 
The values of the dimensionless parameters for the plot were if = 0.072, j3 = 0.104. 
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